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The formation of coherently strained three-dimensional islands on top of the wetting layer in 
Stranski-Krastanov mode of growth is considered in a model in 1+1 dimensions accounting for the 
anharmonicity and non-convexity of the real interatomic forces. It is shown that coherent 3D islands 
can be expected to form in compressed rather than in expanded overlayers beyond a critical lattice 
misfit. In the latter case the classical Stranski-Krastanov growth is expected to occur because the 
misfit dislocations can become energetically favored at smaller island sizes. The thermodynamic 
recison for coherent 3D islanding is the incomplete wetting owing to the weaker adhesion of the edge 
atoms. Monolayer height islands with a critical size appear as necessary precursors of the 3D islands. 
The latter explains the experimentally observed narrow size distribution of the 3D islands. The 2D- 
3D transformation takes place by consecutive rearrangements of mono- to bilayer, bi- to trilayer 
islands, etc., after exceeding the corresponding critical sizes. The rearrangements are initiated by 
nucleation events each next one requiring to overcome a lower energetic barrier. The model is in 
good qualitative agreement with available experimental observations. 

PACS numbers: 68.35.Bs, 68.55.Jk 
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I. INTRODUCTION 

The preparation of arrays of defect free three-dimen- 
sional {3D) nanoscale islands is a subject of intense re- 
search in the last decade owing to possible optoelectronic 
applications as quantum dots. The latter are pronji, 
ing for fabrication of lasers and light emitting diodesEI" 
In recent time the instability of two-dimensional (2D) 
growth against the formation of coherently strained 3D 
islands in highly mismatched heteroepitaxial systems has 
been successfully used to produce quantum dots. This is 
the well known Stranski-Krastanov (SK) growth mode 
where the decrement of the strain energy in the 3D is- 
lands overcompensates the contribution of the surface en- 
ergy. 

When the adhesion forces between the substrate and 
film materials overcompensate the strain energy stored 
in the overlayer owing to the lattice mismatch, a thin 
pseudomorphous wetting layer consisting of an integer 
number of monolayers is first formed by a layer-by-layer 
mode of growth. This kind of growth cannot continue 
indefinitely because of the accumulation of strain energy 
and the disappearance of the energetic influence of the 
substrate after several atomic diameters. Then, in the 
thermodynamic limit, unstrained 3D islands are formed 
and grow on top of the wetting layer, the lattice mis- 
fit being accommodated by misfit dislocations (MDs) at 
the wetting layer - 3D islands boundaryO'cl. Thus the wet- 
ting layer and thCrSD islands represent different phases in 
the sense of Gibbg3 separated by an interphase boundary. 
The energy of the latter is given by the energy of the array 
of MDs. This iS|_the classical Stranski-Krastanov mech- 
anism of growth! (see Fig. ^)). However, it has been 
found that under certain conditions coherently strained 
(dislocation free) 3D islands are formed on top of the 
wetting layer (Fig. (|l|b)). These islands are strained to 
fit the wetting layer in the middle of their bases but, 
more or less strain-free near their top and side wallsQ'li^ 
Such coherently strained islands are formed at large pos- 
itive misfits when the lattice parameter of the overlayer 
is larger than that of the substrate and the overlayer is 
compressed. It has also been observed that the size dis- 
tribution of the 3D islands is very narrow. The above 
observations fe^V^ been reported for the -growth of Ge 
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and InP on Ino.sGao.sPc^ 
lattice misfit is positive and very large (4.2, 7.2, and 
«3.8% for Ge/Si, InAs/GaAs, and InP/Ino.sGao.sP, re- 
spectively) for semiconductor materials which are char- 
acterized by directional and brittle chemical bonds. The 
only exception to the authors' knowledge is the system 
PbSe/PbTe(lll) in which the-piisfit is negative (-5.5%) 
and the,overlayer is expandedEj. However, the authors of 
Ref. (E2I) note that whereas the in-plane lattice parame- 
ter of the PbSc wetting layer is strained to fit exactly the 
PbTe substrate, the parameter of the 3D islands rapidly 
decreases, reaching 95% of the bulk PbSe lattice constant 
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at about 4 monolayers coverag One could speculate 
that the lattice misfit is accommodated by MDs intro- 
duced at the onset of the 3D islanding. 

Whereas the classical SK growth is more or less clear 
from both thermodynamic and kinetic points of view, the 
formation of coherent 3D islands still lacks satisfactory 
explanation. We can consider as a first approximation 
the formation of coherent 3D islands in SK growth as 
homoepitaxial growth on uniformly strained crystal sur- 
face both film and substrate materials having one and the 
same bonding. If so, it is not clear what is the thermo- 
dynamic driving force for 3D islanding if the islands are 
coherently strained to the same degree as the wetting 
layer. It is also not clear why coherent 3D islands are 
observed in compressed rather than in expanded overlay- 
ers. Another question which should be answered is why 
the formation of coherent 3D islands requires very large 
value of the positive misfit. The reason for the narrow 
size distribution is still unclear although much effort has 
been made to elucidate the problemcSEj. Finally, the 
mechanism of formation of coherent 3D islands is still an 
open question. 

Two major approximations are usually made when 
dealing theoretically with the formation of coherently 
strained 3D islands. The first is the use of the linear 
theory of elasticity in order to compute thC|-sti-pn^(m- 
tribution to the total energy of the islandtflMBEl M. 
However, the validity of the latter is hard to accept bear- 
ing in mind the high values of the lattice mismatch. As 
will be shown below the MDs differ drastically in com- 
pressed and expanded films. Second, it is commonly ac- 
cepted that the interfacial energy between the wetting 
layer and the dislocation free 3D islands is sufficiently 
small and can be neglected in the case of coherent SK 
growth. The latter is equivalent to the assumption that 
the substtate-Ithe-wetting layer) wets completely the 3D 
islandsE3'Lll"L3E3"E3. In fact this assumption rules out 
the 3D islanding from thermodynamic point of view as 
3D islands are only possible at incomplete wetting, or in 
othMJSfocds, when the interfacial energy is greater than 
zercOLJLj. As shown below the adhesion of the atoms 
to the wetting layer is also distributed along the island in 
addition to the strain distribution and plays a more sig- 
nificant role than the latter. Due to the lattice misfit the 
atoms are displaced from their equilibrium positions in 
the bottoms of the potential troughs they should ocupy 
at zero misfit. In such a way the adhesion of the atoms 
to the substrate is stronger in the middle of the islands 
and weaker at the free edges. The average adhesion of 
an island of a finite size is thus weaker compared with 
that of an infinite monolayer. An interfacial boundary 
appears and the wetting of the island by the substrate 
(the wetting layer) is incomplete in the average. It is this 
incomplete wetting which drives the formation of dislo- 
cation free 3D islands on the uniformly strained wetting 
layer. 

In the present paper we make use of a more realistic 
interatomic potential which is characterized by its an- 
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harmonicity, in the sense that the repulsive branch is 
steeper than the attractive branch, and by its noncon- 
vexity, which means that it possesses an inflection point 
beyond which its curvature becomes negative. Recently 
Tarn and Lam have used a Mie potential to describe the 
mode of growth in a kinetic Monte Carlo procedureCJ. 
However, the above mentioned authors did not study the 
effect of misfit sign. Moreover, the distribution of the 
stress in the 3D islands has beea studied again withm 
the continuum elasticity theoryc^. Yu and MadhukarO 
computed the energy and the distribution of strain in co- 
herent Ge islands on Si(OOl) using a molecular dynamics 
coupled with the Stillinger- Weber potentialcZl but did not 
study the effect of anharmonicity in the general case. 

The use of such a potential allows us to answer the 
question why coherently strained 3D islands appear pre- 
dominantly in compressed overlayers. Comparing the en- 
ergies of mono- and multilayer islands allows to make 
definite conclusions concerning the mechanism of forma- 
tion and growth of the 3D islands, and the thermody- 
namic reason for the narrow size distribution. It turns 
out that there is a critical 2D island size above which 
the monolayer islands become unstable against the bi- 
layer islands. Thus,|-as has been shown earlier by Stoy- 
anov and MarkovBH, Priester and LannooEa and Chen 
and Washburrfii, the monolayer islands appear as neces- 
sary precursors for the formation of 3D islands. Beyond 
another critical size the bilayer islands become unstable 
against trilayer islands, etc. Then, the growth of 3D is- 
lands consists of consecutive transformations. As a result 
of each one of them the islands thicken by one monolayer. 
The critical size for the mono-bilayer transformations in- 
creases sharply with the decrease of the lattice misfit go- 
ing asymptotically to infinity at some critical misfit. The 
monolayer islands are thus always stable against the mul- 
tilayer islands below this critical misfit which explains the 
necessity of large misfit in order to grow coherent 3D is- 
lands. The critical misfit in expanded overlayers is nearly 
twice greater in absolute value than that in compressed 
overlayers which in turn explains why coherent 3D is- 
landing is very rarely (if at all) observed in expanded 
overlayers. 

The edge atoms are more weakly bound than the atoms 
in the middle of the islands. This is due to the weaker 
adhesion of the edge atoms to the wetting layer. Thus, 
the 2D-3D transformation takes place by transport of 
atoms from the edges of the monolayer islands, where 
they are weakly bound, on top of their surfaces to form 
islands, p£ the upper layer where they are more strongly 
bouncl3'C3. This process is then repeated in the transfor- 
mation of bilayer to trilayer islands, etc. The critical size 
for the 2D-3D transformation to occur is the thermody- 
namic reason for the narrow size distribution of the 3D 
islands. 

In the case of expanded overlayers the atoms interact 
with each other through the weaker attractive branch of 
the potential and most of the atoms are not displaced 
from their equilibrium positions. The size effect is very 



weak, the average adhesion is sufficiently strong, and the 
critical sizes for 2D-3D transformation either do not exist 
or appear under extreme conditions of very large abso- 
lute value of the misfit. In any case MDs are introduced 
before the formation of bilayer islands. The coherent 
monolayer islands are either energetically stable against 
multilayer islands or MDs are introduced before the 2D- 
3D transformation. As a result the classical SK growth 
is expected in expanded overlayers. 



II. MODEL 

We consider a model in 1+1 dimensions (substrate -I- 
height) which we treat as a cross-section of the real 2-|-l 
case. An implicit assumption is that in the real 2-1-1 
case the monolayer islands have a compact rather than a 
fractal shape and the lattice misfit is one and the same 
in both orthogonal directions. Although the model is 
qualitative it gives correctly all the essential properties 
of the reaL-2-tl system as shown by Snyman and van 
der MerwecjEil. In this model the monolayer island is 
represented by a finite discrete Frenkel-Kontorova linear 
chain of atoms subject to an external-^eHodic potential 
exerted by a rigid substrate (Fig. (P)E3~Eil. We consider 
as a substrate the uniformly strained wetting layer of the 
same material consisting of an integer number of mono- 
layers. In other words, we consider the SK growth in two 
separate stages. The first stage is a Frank-van der Merwe 
(layer-by-layer) growth during which the wetting layer is 
formed. The second stage is a Volmer- Weber growth of 
3D islands on top of the wetting layer. In this paper we 
restrict ourselves to the consideration of the second stage 
assuming the wetting layer is already built up. The ener- 
getic influence of the initial substrate is already lost and 
the bonding between the atoms in the 3D islands is the 
same as that of the atoms of the first atomic plane of 
the 3D islands to the atoms belonging to the uppermost 
plane of the wetting layer. 

The atoms of the chain are connected rSSiith bonds that 
obey the generalized Morse potentialP^'" 
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shown in Fig. (Q) where ^ and v {^j. > v) are constants 
that govern the repulsive and the attractive branches, 
respectively, and b is the equilibrium atom separation. 
For ^ = 2v the potential (^ turns into the familiar Morse 
potential. In the case of homoepitaxy the bond strength 
Vo is related to the energy barrier for desorption. 

The potential (^ possesses an infiection point Xmf = 
b + ln(/i/j/)/(/i — v) beyond which its curvature becomes 
negative. The latter leads to a distortion of the inter- 
atomic bonds in the sense that long, weak and short, 
strong bonds, alternate (see the upper right-hand corner 
of Fig. (0)EjEa), and to the appearance of structures 
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consisting of nuJtiple MDs (multikinks or kink-antikink- 
kink solutions)LZI. The latter represent two kinks (or soli- 
tons) connected by a strongly stretched out bond (the 
antikink) . 

The 3D islands can be represented by linear chains 
stacked one upon the other as-jp the model proposed 
by Stoop, and van der Merwa^, and by Ratsch and 
Zangwilo, each upper chain being shorter than the lower 
one. In principle, the Frenkel-Kontorova model is inad- 
equate to describe a thickening overlayer because of two 
basic assumptions inherent in it. The first one is the 
rigidity of the substrate. Assuming that the substrate 
remains rigid upon formation of 3D islands on top of 
it rules out the interaction between the islands through 
the elastic fields around them. It is believed that this 
assumption is valid for very thin deposits not exceeding 
one or two monolayers. The second one is connected with 
the relaxation effects. When a new monolayer island is 
formed on top of the previous one the latter should re- 
lax and the strains in the island will redistribute. One 
can expext that the formation, say, of a second mono- 
layer will make the bonds between the first monolayer 
atoms effectively stiffer. As will be discussed below this 
will lead to weaker adherence of the atoms in the first 
monolayer to the wetting layer. MDs could be also in- 
troduced to relieve the strain. Nevertheless, the Frenkel- 
Kontorova model can provide excellent qualitalrtie-gen- 
eralizatiop,in two dimensions both horizontally&O, and 
verticallyEJ. According to the authors of Ref. (e3) n- 
layer island can be mimicked by assuming that the force 
constant of the interatomic bonds is n times greater than 
that of a monolayer island. Thus a bilayer island under 
compression could be simulated by doubling the value of 
the repulsive constant fi. This approach obviously gives 
the upper bound of the effect of the next layers on the 
redistribution of the strain in the lower layers. An im- 
plicit shortcoming of this method is that it assumes the 
same number of bonds (and correspondingly atoms) in 
the upper chains and thus does not allow calculations of 
clusters with different slopes of the side walls. 

Another approach tq-the problem has been proposed 
by Ratsch and ZangwillEJ. They accepted that each layer 
(chain) presents a rigid sinusoidal potential to the chain 
of atoms on top of it. The atom, or more precisely, the 
potential trough separation of the lower chain is taken 
as average of all atom separations. As the strains of the 
bonds which are closer to the free ends are smaller, the 
average atom separation 6„ in the nth chain is closer 
to the unperturbed atom spacing b and the lattice misfit 
fn+i = {b—bn)/bn for computing the energy of the n+lst 
chain is smaller in absolute value than the misfit f = {b~ 
a) I a which is valid only for the base chain that is formed 
on the wetting layer the latter having an atom separation 
a. In such a way the lattice misfit and in turn the bond 
strains gradually decrease with the islands height. Every 
upper chain was taken shorter than the lower one by an 
arbitrary number of atoms and was centered on top of 
it as shown schematically in Fig. (|3|). Moreover, every 



uppermost chain is taken frozen (the relaxation of the 
lower chain upon formation of the next one is ruled out) 
and serves as a template for the formation of the next 
one. Then, the formation of each next chain does not 
exert any influence on the distribution of strain in the 
previous chains, and thus this approach represents the 
lower bound of the effect of the next layer. 

In the preseat paper we will use the approach of Ratsch 
and Zangwillc3. The main reason is that it allows a grad- 
ual attenuation of the strain with the island height, and 
also different angles of the side walls. We believe that 
although rather crude this approach gives correctly the 
essential physics with one exception. It does not account 
for the decrease of the average adhesion of the base chain 
to the wetting layer upon thickening of the islands. An 
approximate evaluation of the latter effect can be ob- 
tained by using the upper bound approach. It should 
be emphasized that both approaches show qualitatively 
identical vesnlts. We could expect that the results of more 
accurate calculations including the strain relaxation will 
not differ qualitatively by those presented below. Pre- 
liminary studies with an energy minimization program 
allowing strain relaxation always produced dislocated ex- 
panded and coherent compressed islands in agreement 
with the results shown below. Note that owing to the 
approximations of the model (1-1-1 dimensions and the 
lack of relaxation) the figures obtained as a result of the 
calculations, e.g. 3.25% for the critical misfit for 3D is- 
landing, should not be taken as meaningful. Finally, we 
have to mention that the numerical solution of the sys- 
tem of governing equations (^ requires no more than 
a seconds on a 100 MHz PC even when the number of 
equations (atoms in the chain) is about 100. 

Discussing the stability of mono- and multilayer is- 
lands wp-foUow the approach developed by Stoyanov and 
Markovca. We start from the classical concept of mini- 
mum of rtte surface energy at a fixed volume. Following 
StranskiEil, the surface energy F{N) is defined as the dif- 
ference between the potential energy of the cluster con- 
sisting of N atoms and the potential energy of the same 
number of atoms in the bulk crystal 

N 

F{N) Nct>k - '^•^ 

i=l 

which is valid for clusters with arbitrary shape and size. 
Here 0^ is the work necessary to detach one atom from 
a kink position (or the energy of atom in the bulk of 
the crystal), and the sum gives the work required to dis- 
integrate the cluster into single atoms. Since the term 
iV(/)fc does not depend on the cluster shape the stabil- 
ity of mono- and multilayer islands is determined by the 
above sum. The maximum of the sum corresponds to 
a minimum of the surface (edge) energy of the cluster. 
Therefore, as a measure of stability, we adopt the po- 
tential energy per atom of the clusters, which is, in fact, 
equal to the above sum taken with a negative sign. 
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The potential energy of a chain of the nth layer con- 
sisting of Nn atoms reads 



Af„-1 



En 



where 
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(3) 



accounts for the adhesion of the ith atom. Xi are the 
coordinates of the atoms taken from an arbitrary origin. 
The difference AXi — X^+i — Xi is in fact the distance 
between the i + 1st and zth atoms. The first sum in Eq. 

gives the energy of the bond strains. The second sum 
gives the energy of the atoms in the periodic potential 
field created by the lower chain where W is its amplitude 
and bn-i is the average potential trough separation of 
the underlying layer. In general W should be a function 
of the atom separation of the underlying layer and thus 
should depend on n but for simplicity we neglect this 
dependence. As mentioned above = a holds only 

for the base chain. The amplitude W can be considered 
in our model as the barrier for surface diffusion. On a 
nearest neighbor bond hypothesis W is related to the 
substrate-deposit bond strength by 



W = gVo 



(4) 



where <? < 1 is a constant of proportionality varying 
approximately from 1/30 for long-range vaxL-der Waals 
forces to 1 /3 for short-range covalent bondso. 

The average of the second sum in Eq. (^) for the base 
chain divided by Vo 
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has the same physical meaning as the adhesion parameter 
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2a 
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which accounts for the incomplete wetting of the 3D is- 
lands by the substrate in heteroepitaxy (tr, ai and as 
being the specific surface energies of the overlayer, the 
interface and the substraije, respectively, and p being the 
specific adhesion energy )eI. In the case of the classical SK 
growth the adhesion parameter is given bjt, $ = ed/2a 
where is the energy of a net of MDsEj. We have 
the case of complete wetting when < 0. The forma- 
tion of 3D, islands can obviously take place only when 

< <z> < ila. 

Minimization of En with respect to Xi results in a set 
of governing equations for the atom coordinates in the 
form 



Asin(27r6) = 0, (6) 



where = 6„_i(^i — — /„) is the strain of the 
ith bond, = Xi/bn^i is the displacement of the zth 
atom with respect to the bottom of the ith potential 
trough, /„ is the misfit between the nth chain and the 
substrate potential exerted by the n — 1st chain, and 
A = ttW{ii — v) / fjLvhn-iVo- The lattice misfit has its 
largest value / = {b — a)/a only for the base chain in mul- 
tilayer islands, and goes to zero with increasing islands 
thickness. Expanding of the exponentials in Taylor series 
for small strains gives the se t . of . equations that govern 
the discrete harmonic modelr^ri. Solving numerically 
the system of equations gives the atom displacements 
^i and all the parameters characterizing the system can 
be easily computed. 

The properties of the solutions of the system are 
of crucial importance for understanding of the coherent 
SK growth. Two forces act on each atom: first this is 
the force exerted by the neighboring atoms, and second, 
this is the force exerted by the substrate (the underlying 
chain or the wetting layer). The first force tends to pre- 
serve the natural spacing b between the atoms, whereas 
the second force tends to place all the atoms at the bot- 
toms of the corresponding potential troughs of the sub- 
strate separated at a distance b. As a result of the 
competition between the two forces the bond strains and 
the atom adhesion are distributed along the chain. The 
undislocated solution (Fig. (|^a)) clearly shows the de- 
crease of the atom adhesion at the ends of the chain as 
the atoms are more and more displaced towards the chain 
ends. In the case of positive misfit the dislocation repre- 
sents an empty potential trough the bond in the core of 
the dislocation being strongly stretched out (Fig. (^)). 
This picture is equivalent to a crystal plane in excess in 
the substrate. In the opposite case of negative misfit (Fig. 
(0c)) the dislocation represents two atoms in one trough 
(a crystal plane in excess in the overlayer), the bond in 
the dislocation core being compressed. Both configura- 
tions are energetically equivalent in the harmonic approx- 
imation where the force between the atoms increases lin- 
early with the atom separation. This is not, however, 
the case when an anharmonic potential is adopted. The 
latter displays a maximum force between the atoms at 
X = Xinf- This is the theoretical tensile stress of the 
material atens = Vajjiiy / and if the actual force 
exerted on the correspoBjdiag-iond is greater than atens 
the bond will break upEa^EZICj. Thus the interval of ex- 
istence of dislocated solutions in compressed chains de- 
pends on the material parameters Vo, W, fj,, v, /, and be- 
comes very narrow. Dislocated solutions ia-compressed 
chains exist only in sufficiently long chainsE^tJ beyond 
some critical chain length. As will be shown below, this 
leads to coherent SK growth in compressed overlayers. 
On the contrary, the bonds in the cores of the MDs in 
expanded chains are compressed and cannot break. As 
a result MDs become energetically favored and can be 
introduced in very short chains. Thus, the classical SK 
growth should be expected in expanded overlayers as the 



5 



dislocated islands with a monolayer height can become 
energetically favored long before the coherently strained 
multilayer islands. 

III. RESULTS 
A. Monolayer islands 

The distribution of the bond strains along the chains 
is shown in Fig. (|^a). As expected the bonds in the 
middle of the chains are strained to fit exactly the uni- 
formly strained wetting layer. The strains at the chain 
ends tend to zero. In fact the strains of the hypothetiCip.- 
roth and A^th bonds should be exactly equal to zeroHo. 
The strains in the middle of the expanded chain com- 
pared with these of compressed ones are much closer to 
— / owing to the weaker attraction between the atoms 
of the chain. Fig. (^) shows the distribution of the 
bond energy. It is seen that in the case of compressed 
chains (/ > 0) the bond energy in the middle of the 
chain is smaller than that in expanded chains owing to 
the stronger atom repulsion. 

The distribution of the adhesion of the separate atoms 
<l>i (Eq. (^)) (taken in terms of the bond energy Vo as 
(fPi ~ Vo)/Vo) is demonstrated in Fig. (^). The weaker 
adhesion at the chain ends, which is often overlooked 
in theoretical models, is due to the displacement of the 
atoms from the bottoms of the potential troughs (see 
Fig. (^a)). What is more important is that the atoms in 
the expanded chains adhere much more strongly to the 
wetting layer compared with the atoms in the compressed 
chains. 

Fig. (0) shows the dependence of the mean adhesion 
parameter <P (Eq. (||)) on the number of atoms. As can 
be expected the atom adhesion in expanded overlayers 
is stronger than that in compressed ones owing to the 
weaker attraction between the atoms in the former. The 
forces exerted from the substrate are stronger than the 
forces between the chain atoms and the latter are situated 
more deeply in the potential troughs. The curves display 
maxima which are due to the interplay between the frac- 
tion of the most strongly displaced end atoms and the 
values of the particular displacements. In short chains 
the atoms are weakly displaced from the bottoms of the 
potential troughs and the adhesion is stronger. With in- 
creasing chain length the displacements of the end atoms 
increase and beyond some length saturate and do not in- 
crease anymore. The fraction of weakly displaced middle 
atoms increases and a maximum is displayed. The value 
of the maximum (not shown) decreases sharply with de- 
creasing misfit going asymptotically to zero at zero misfit. 
This means that ^ > at any value of the nonzero misfit 
which is the thermodynamic reason for 3D islanding. 

Fig. (^) shows the distribution of the total energy 
(strain plus adhesion) in chains with positive and neg- 
ative misfit. The maxima in the middle are due to the 



strain contribution whereas the increase of the energy at 
the ends is due to the weaker adhesion. It is first seen that 
the atoms in the expanded chain are considerably more 
strongly bound to each other and to the substrate. The 
main difference between both curves is that the atoms at 
the free ends in compressed chains are much more weakly 
bound than the end atoms in expanded chains. This re- 
sult is of crucial importance for our understanding of the 
mechanism of transformation of the mono- to multilayer 
(3D) islands. We conclude from Fig. ^ that compressed 
islands display a greater tendency to transform into bi- 
layer islands and further to form coherent 3D islands in 
comparison with expanded islands. 

B. Multilayer islands 

The multilayer (3D) islands can be full or frusta of 
pyramids and can have side walls with different slope. 
The effect of the side walls slope on the minimum en- 
ergy shape is more or less clear. More unsaturated dan- 
gling bonds normal to the film plane appear on side walls 
with smaller slope and the corresponding surface energy 
is greater. Obviously, the surface energy of the steepest 
walls with a slope of 60° is the lowest one. One could 
expect that the islands bounded with the steepest walls 
will be more stable than the flatter islands. The problem 
whether the pyramids are full or frusta is more difficult 
to resolve. First, with increasing pyramid height the lat- 
tice misfit decreases and the mean strain vanishes. This 
in turn leads to increase of the adhesion of the separate 
atoms, and, as a whole, to an increase of the bond en- 
ergy closer to the apex of the pyramids. On the other 
hand, the layers which are closer to the apex are smaller 
in size and the size effect increases. The latter leads to 
smaller work of evaporation per atom of a whole upper- 
most atomic plane. As has been known for a long time, 
the work required to disintegrate a whole atomic plane 
into single atoms (the mean separation work) taken with 
a negative sigHpis equal to its chemical potential at the 
absolute zeroOO. Hence, adding to the pyramid smaller 
and smaller upper base atomic planes leads to a decrease 
of the mean separation work of the upper base and in 
turn to higher chemical potential. As a result we could 
expect that frusta of pyramids with a slope of 60° of the 
side walls will be energetically favored. This is clearly 
seen in Fig. (^ which demonstrates the energy per atom 
of pyramids with different side walls slopes as a func- 
tion of the height taken as a number of monolayers. The 
curves display minima at a certain height which clearly 
show that the frusta of the pyramids are the lowest energy 
configurations. The energy of the full pyramids is much 
higher. The minimum of the 60° side wall slope is the 
lowest one thus confirming the above consideration. The 
steepest side wall slope of 60° is a natural consequence 
of the model which considers a face centered cubic rather 
than a diamond lattice. It is worth noting that Ratsch 
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and Zanguill also report that the steepest side walls are 
energetically favoredcfl. 

The above result does not mean that in the real ex- 
periment the coherent 3D islands will grow as frusta of 
pyramids. The lowest minimum in Fig. (^) represents in 
fact the equilibrium shape of the islands. In reality the 
crystallites grow with a shape which is determined by the 
rates of growth of the-|different walls and thus depends on 
the supersaturatiorO. The growing crystal is bounded 
by the walls with the lowest growth rate at the given su- 
persaturation. Mo et al have established with the help of 
scanning tunneling microscopy (STM) that small coher- 
ently strained Ge islands ("hut" islands) grow-on Si(OOl) 
as full pyramids bounded with (105) side wallsllil, whereas 
Voigtlander and Zinner observed frusta of tetrahedron 
Ge pyramids on Si(lll) with aspect (height-to-base) ra- 
tio showing a maximum of about 0.135 at a coverage of 
4 MLstj. All the above is valid for sufficiently large crys- 
tals. We are interested here of the initial stages of growth 
of the 3D islands, or more precisely, in the transforma- 
tion of monolayer into multilayer islands. As shown in 
the next sections the formation and growth of 3D islands 
proceeds by consecutive transformations of monolayer is- 
lands into bilayer and then into multilayer islands which 
is the lowest energy path of the 2D-3D transformation. 

It should be stressed that the adhesion parameter of 
a monolayer island should differ significantly from that of 
a multilayer island with the same base chain length. In 
our model they are equal. The reason is that the model 
does not allow the relaxation of the lower chains after 
formation of new ones on top of them. The above is obvi- 
ously incorrect as the formation of a second chain on top 
of the base one leads_tp effectively stronger lateral bonds 
in the bilayer islandsEj. We will try to qualitatively eval- 
uate this problem and to discuss its consequences. As 
mentioned above the bilayer island could be treated as 
a first approximation as a monolayer island with a dou- 
bled force constantESl. As a result both the fraction of the 
strongly displaced end atoms and the corresponding dis- 
placements will be larger. Then the adhesion parameter 
of a bilayer island will be greater than that of a monolayer 
island with the same width. An evaluation of this effect 
can be made by using the approach of van der Merwe et al 
mentioned aboveE3 by doubling of the constant ^ in com- 
pressed chains. Thus for mono-, bi- and trilayer islands 
with ^ = 12, 24 and 36, one obtains 0.024, 0.066 

and 0.1, respectively (z^ = 6, / — 0.05, N — 21). As seen 
the effect of the third layer is weaker than that of the 
second which is easy to understand. The effect of for- 
mation of the next monolayers will have a smaller effect 
on the adhesion of the island and after some thickness 
the adhesion parameter will not change anymore. Thus 
the base layer atoms in a coherent multilayer island are 
more weakly bound to the wetting layer. What follows is 
that once formed the bilayer islands stabilize the further 
growth of the coherent 3D islands. 



C. Stability of mono- and multilayer islands 

We compare further the energies of mono- and mul- 
tilayer islands with different thickness. The latter are 
bounded with 60° side walls as they have the lowest min- 
imum energy as shown above. Fig. ([lO|a) shows the 
dependence of the energy of compressed monolayer and 
multilayer islands on the total number of atoms at com- 
paratively small lattice misfit of 3%. As seen the mono- 
layer islands are always stable against bilayer and trilayer 
islands. A 2D-3D transformation is thus not expected 
and the film should continue to grow in a layer-by-layer 
mode coupled with an introduction of MDs at a later 
stage. The same dependence but at a larger misfit of 
5% is demonstrated in Fig. (p^). The monolayer is- 
lands become unstable against the bilayer islands beyond 
a critical island size N12, the bilayer islands in turn be- 
come unstable against the trilayer islands beyond a sec- 
ond critical number A^23, etc. The curve denoted by MD 
represents the energy of a monolayer chain containing 
one MD. The latter begins at a large number of atoms 
{N = 52) because the bonds in the cores of the MDs 
break up for shorter chains. This is due to the fact that 
the force exerted on these bonds from the neighboring 
atoms is greater than the theoretical tensile stress of the 
film material atens as mentioned above. Curve 1 which 
represents the energy of undislocated monolayer chain is 
computed for clarity up to a number of atoms smaller 
than the number (52) at which the solutions of the dislo- 
cated chain appear. The reason is that the values of the 
energy are very close and the curves are undistinguish- 
able for the eye. The energies of monolayer chains with 
and without MDs cross each other at about N = 300 
(not shown) which means that coherent 3D islands are 
formed long before the introduction of MDs. Moreover, 
the dislocated chain with a monolayer height has an en- 
ergy much higher than the energies of the undislocated 
multilayer islands. The latter clearly shows that the film 
"prefers" to grow as coherent 3D islands in which the 
gradual decrease of the strain energy overcompcnsates 
the surface energy rather than to introduce MDs in the 
first monolayer. 

Fig. ( [TT| ) demonstrates the same dependence as in Fig. 
(|l0| ) but in expanded chains. The absolute value of the 
negative misfit is very large (-10%). At absolute values 
of the misfit smaller than 5.5% (not shown) the behavior 
of the energies is the same as in Fig. (^0|a). The energies 
of the coherent mono- and multilayer chains cross again 
each other at some critical numbers of atoms but the 
dislocated monolayer chain (denoted by MD) becomes 
energetically favoured noticeably before the coherent bi- 
layer chain to become stable. The classical SK growth 
should take place in expanded overlayers. 

Fig. dl^ ) shows the misfit dependence of the first 
critical size N12 for both positive and negative misfits. 
As seen it increases sharply with decreasing misfit going 
asymptotically to infinity at some critical misfits denoted 
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by the vertical dashed hnes. The existence of a critical 
positive misfit for coherent SK growth to occur explains 
why a high mismatch epitaxy is required in order to grow 
coherent 3D islands. The critical misfit below which the 
expanded monolayer islands arc always stable against 
multilayer islands is nearly twice larger in absolute value 
compared with the same quantity in compressed overlay- 
ers. Thus coherent SK growth in expanded overlayers 
could be observed at unrealistically large absolute values 
of the negative misfit. 

We conclude that the classical SK growth or a 2D 
growth will be observed in the thermodynamic limit at 
small positive misfits and a coherent SK growth at misfits 
greater than a critical misfit. This result clearly explains 
why large positive misfit is required for the coherent SK 
growth to occur. The large positive misfit leads to large 
atom displacements and in turn to weaker adhesion. The 
physics is essentially the same as in the case of heteroepi- 
taxial growth of 3D islands directly on top of thp-surface 
of the foreign substrate (Volmcr- Weber growth)c3. 

D. Mechanism of 2D-3D transformation 



is in equilibrium is equal to the difference of the energies 
of desorption of the atoms from the same and the foreign 
substrate, This is namely the driving force for the 2D- 
3D transformation to occur. The 1+1 model is in fact 
one-dimensional and the nuclei do not exist in the ther- 
modynamic sense because the length of a. mw of atoms 
does not depend on the supersaturationEJ'Ell. However, 
considering our 1-1-1 model as a cross-section of the real 
2-1-1 case we can treat the curve 1-2 in Fig. ( p^ ) as the 
size dependence of the free energy for nucleus formation 
and growth. We would like to emphasize that in the 
2+1 case the nucleus does not necessarily consist of one 
atom. Its size should depend on the lattice misfit, and 
in the real situation - on the temperature. The curves 
denoted by 2-3 and 3-4 in Fig. ( [l3|) represent the energy 
changes of bilayer to trilayer islands, and of trilayer to 
fourlayer islands, respectively. As seen they behave in 
the same way and the work for nucleus formation (the 
respective maxima) decrease with the thickening of the 
islands. The latter means that the mono-bilayer trans- 
formation is the rate determining process of the total 
mono-multilayer (2D-3D) transformation. 



It is natural to assume that once the monolayer islands 
become unstable against the bilayer islands {N > N12) 
the former should rearrange themselves into bilayer is- 
lands. As shown below the mono-bilayer transformation 
can be considered as the first step for building sufhciently 
high 3D crystallites. The mechanism of the mono-bilayer 
transformation is easy to predict having in mind that the 
edge atoms are more weakly bound than the atoms in the 
middle. The edge atoms can detach and difuse on top of 
the monolayer islands giving rise of clusters of the second 
layer. We consider first in more detail the transformation 
of a monolayer island (chain) with a length No > N12 
into bilayer island. For this aim we plot the energy E{n) 
of an incomplete bilayer island which consists of No — n 
atoms in the lower layer and n atoms in the upper layer 
referred to the energy Eo of the initial chain consisting 
of No atoms as a function of the number of atoms n in 
the upper layer. This is the curve denoted by 1-2 in Fig. 
(|l3|). As seen it displays a maximum at n = 1 after 
which AEn = E{n) — Eo decreases upto the complete 
mono-bilayer transformation at which n = {No — l)/2. 

Curve 1-2 in Fig. (^^ has the characteristic behavior 
of a nucleation process. The cluster at which the maxi- 
mum of Ai? is observed can be considered as the |esitical 
nucleus of the second layer. As shown in Ref. (e3) the 
mono-bilayer transformation is a real nucleation process 
when 2+1 heteroepitaxial Volmer- Weber model is con- 
sidered, in other words, when the 3D islands arc formed 
directly on top of the foreign substrate without the for- 
mation of an intermediate wetting layer. The chemical 
potential of the upper island at the maximum is exactly 
equal to that of the initial monolayer island, and the su- 
persaturation with which the nucleus of the second layer 



IV. DISCUSSION 

The Stranski-Krastanov growth mode appears as a re- 
sult of the interplay of the film-substrate bonding, strain 
and surface energies. A wetting layer is first formed on 
top of which 3D islands nucleate and grow. The 3D is- 
lands and the wetting layer represent necessarily different 
phases. If this was not the case the growth should con- 
tinue by 2D layers. Then we can consider as a useful 
approximation the 3D islanding on top of the uniformly 
strained wetting layer as Volmer- Weber growth. The lat- 
ter requires the adhesion of the atoms to the substrate 
to be smaller than the cohesion between the overlayer 
atoms. In other words, the wetting of the substrate by 
the overlayer should be incomplete. In the classical SK 
growth this condition is fulfilled because of the forma- 
tion of an array of misfit dislocations at the boundary 
between the islands and the wetting layer. The atoms 
are displaced from the bottoms of the potential troughs 
(mostly in the cores of the MDs, see Fig. (|b)) and (|c)) 
and thus are more weakly bound in average to the un- 
derlying wetting layer, irrespective of that the chemical 
bonding is one and the same. As a result the lattice 
misfit gives rise to an effective adhesion which is weaker 
than the cohesion of the overlayer atoms. Contrary to 
the wetting layer the 3D islands are elastically relaxed 
and their atom density differs from that of the former. 
Thus, the wetting layer and the 3D islands really rep- 
resent different phases separated by a clear interfacial 
boundary whose energy is in fact the energy of the ar- 
ray of MDs. The physical reason for 3D islanding in the 
coherent SK growth is practically the same. In this case 
the atoms near the islands edges are displaced from the 
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bottoms of the corresponding potential troughs (see Fig. 
(^)) and they adhere more weakly to the wetting layer 
compared with the atoms in the middle. The thicker the 
islands the stronger is this tendency. Thus, the average 
adhesion of the 3D islands to the wetting layer is again 
weaker than the cohesion in the islands themselves. Then 
we can treat the coherent SK growth as a Volmer- Weber 
growth on top of the wetting layer. The main difference is 
that in Volmer- Weber growth the adhesion parameter <l> 
is constant whereas in the coherent SK mode it depends 
on the islands thickness. 

The weaker adhesion means in fact an incomplete wet- 
ting which appears as the thermodynamic driving force 
for the 3D islanding. The smaller the misfit the smaller 
are the displacements of the edge atoms and in turn the 
stronger is the average wetting. The latter leads to the 
appearance of a critical misfit below which the edge ef- 
fects do not play a significant role. The average wetting is 
very strong and the formation of coherent 3D islands be- 
comes thermodynamically unfavored. The film will con- 
tinue to grow in a 2D mode untill the strain is relaxed by 
introduction of MDs or dislocated 3D islands at a later 
stage. The existence of a critical misfit for 2D-3D trans- 
formation to occur both in compressed and expanded 
overlaj/ers has been noticed in several studies. Pinczolits 
et a£ZI have found that deposition of PbSei_a;Tea; on 
PbTc(lll) remains purely two dimensional when the mis- 
fit is less tharul.6% in absolute value (Se content < 30%). 
Leonard et aa have successfully grown quantum dots of 
In^Gai_^As on GaAs(OOl) with x = 0.5 (/ « 3.6%) but 
60l thick 2D quantum wells at a; = 0.17 (/ w 1.2%). A 
critical misfit of 1.4% has been found by Xie et al upon 
deposition of Sio.sGco.s films on relaxed buffer layers of 
Sij^Gei-a; with varying compositionEa. 

The average adhesion (the wetting) depends strongly 
on the anharmonicity of the interatomic forces. Ex- 
panded islands adhere more strongly to the wetting layer 
and the critical misfit beyond which coherent 3D island- 
ing is possible is much greater in absolute value com- 
pared with that in compressed overlayers. As a result 
coherent SK growth in expanded films could be expected 
at very (unrealistically) large absolute values of the neg- 
ative misfit. The latter, however, depends on the ma- 
terials parameters (degree of anharmonicity, strength of 
the chemical bonds, etc.) of the particulac, system and 
cannot be completely ruled out. Xie et a£3 studied the 
deposition of Sio.sGeo.s films in the whole range of 2% 
tensile misfit to 2% compressive misfit on relaxed buffer 
layers of Si^^Gci-a; starting from x — Q (pure Ge) to a; = 1 
(pure Si) and found that 3D islands are formed only un- 
der compressive misfit larger than 1.4%. Films under 
tensile strain were thus stable against 3D islanding in 
excellent agreement with the predictions of our model. 

The weaker average adhesion in compressed overlayers 
leads to another effect at misfits greater than the critical 
one. At some critical number of atoms N12 the mono- 
layer islands become unstable against the bilayer islands. 
The latter become in turn unstable against trilayer is- 



lands beyond another critical number iV23, and so on. 
As a result the complete 2D-3D transformation should 
take place during growth by consecutive transformations 
of mono- to bilayer, bi- to trilayer islands, etc. Owing to 
the stronger interatomic repulsive forces the edge atoms 
in the compressed monolayer islands adhere more weakly 
to the wetting layer compared with the edge atoms in ex- 
panded islands. This results in an easier transformation 
of mono- to bilayer islands which is the first step to the 
complete 2D-3D transformation. The latter includes also 
kinetics in the sense that the edge atoms have to detach 
and form the upper layers. However, it is not the strain 
at the edges (which is nearly zero) that is responsible 
for the easier detachment of the edge atoms as suggested 
by Kandel and KaxirasE3 but the weaker adhesion. The 
2D-3D transformation is hindered in expanded islands 
as the edge atoms adhere more strongly to the wetting 
layer. On the other hand, the existence of such critical 
sizes, which determine the intervals of stability of islands 
with different thickness, could be considered as the ther- 
modynamic reason for the narrow size distribution of the 
3D islands which is observed in the experiment. The 
latter does not mean that this is the only reason. Elas- 
tic interactions between islands and growth kinetics can 
have greater effect than the thermodynamics. The 2D- 
3D transformation takes place by consecutive nucleation 
events, each next one requiring to overcome a lower en- 
ergetic barrier. Thus, the mono-bilayer transformation 
appears as the rate determining process. 

Let us consider all the above from another point of 
view. The results displayed in Fig. (p^) show that 
the equilibrium shape aspect ratio increases gradually 
with the islands volume. The consecutive stability of 
islands with increasing thickness reflects the fact that 
the increase of the pyramid height is discreet (layer af- 
ter layer) whereas the base chain length remains nearly 
constant. The stronger the adhesion or the smaller the 
misfit the wider will be the intervals of stability of is- 
lands with a fixed height and vice versa. The formation 
of every new crystal plane on the upper crystal face re- 
quires the appearance of a 2D nucleus. As the growing 
surface is usually very small the formation of one nucleus 
is sufficient for the growth of a new crystal plane. Thus 
we could pyppr-t a mononucleus layer-by-layer growth of 
the pyramidsoE3. The latter has been independently es- 
tablished by using of-a kinetic Monte Carlo method by 
Khor and Das Sarmat2l. It should be noted that Duport, 
Priester and Villain established that the monolayer is- 
lands are thermodynamically favored upto a critical size 
beyond whieh, the equilibrium shape becomes nearly a 
full pyramidcZI. The transition from a monolayer island 
to a pyramid is of first order and requires the overcoming 
of an activation barrier which is proportional to /~^. 

It should be stressed that our definition of the critical 
2D island size A''i2 for 2D-3D transformation to besiii 
differs from that in the papjeijs of Priester and Lannoocfl, 
and of Chen and WashburnEll. The former authors define 
the critical size by comparing the energy per atom of 
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monolayer islands with that of fully built 3D pyramids. 
Chen and Washburn have accepted as critical the size 
at which jthe energy of the monolayer islands displays a 
minimumEj. They found also that the critical size Nc 
determined by the minimum of the energy increases very 
steeply with decreasing misfit {Nc oc /^^). Although our 
definition of Nc is different we also observe a very steep 
misfit dependence (see Fig. ([l^)). 

A rearrangement of monolayer height (2D) islands into 
nmjtilayer (3D) islands has been reported by Moison et 
atB who established that the InAs 3D islands begin to 
form on GaAs at a coverage of about 1.75 ML but then 
the latter suddenly decreases to 1.2 ML. This decrease 
of the coverage in the second monolayer could be inter- 
preted as a rearrangement of an amount of nearly half a 
monolayer into 3D islands. The same phenomenon has 
been noticed by Siiklyaev, Shibata and Ichikawa in the 
case of Ge/Si(lll)0. Voigtlander and Zinner noted that 
Ge 3D islands in Ge/Si(lll) epitaxy have been observed 
at the same locations where 2D islands locally jexceeded 
the critical wetting layer thickness of 2 bilayersllj. 

Contrary to the linear theory of elasticity the anhar- 
monicity and the non-convexity of the real interatomic 
potentials lead to different intervals of existence of mis- 
fit dislocations in compressed and expanded overlayers. 
The nonconvexity of the interatomic potential gives rise 
to the possibility of breaking of the expanded bonds in 
the cores of the MDs in compressed overlayers when the 
force exerted on them is greater than the theoretical ten- 
sile strength of the material. As a result MDs in com- 
pressed overlayers appear in sufficiently large islands and 
small coherent 3D islands can appear before that. On 
the contrary, this restriction does not exist in expanded 
overlayers where the bonds in the cores of the MDs are 
compressed. The introduction of MDs can thus become 
energetically favored in short chains (small islands) be- 
fore the formation of coherent 3D islands and the classical 
SK growth should be observed in most cases. 

It should be noted that the results presented above 
depend on the approximations of the model particularly 
when the energy of the multilayer islands is computed. 
Allowing a strain relaxation of lower layers when new 
layers are formed on top of them could lead to earlier 
introduction of MDs but also to weaker adhesion of the 
3D islands to the wetting layer. Thus, applying a more 
refined approach which accounts for the strain relaxation 
in the islands, as well as in the wetting layer, will allow us 
to study the transition from the coherent to the classical 
(dislocated) Stranski-Krastanov growth mode. 

In summary, accounting for the anharmonicity and 
the non-convexity of the real interatomic potentials in 
a model in 1-1-1 dimensions, we have shown that coher- 
ent 3D islands can be formed on the wetting layer in the 
SK mode predominantly in compressed overlayers at suf- 
ficiently large values of the misfit. Coherent 3D islanding 
in expanded overlayers could be expected as an exception 
rather than as a rule. Monolayer height islands with a 
critical size appear as necessary precursors of the 3D is- 



lands. The latter explains the narrow size distribution of 
the 3D islands from thermodynamic point of view. 
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FIG. 1. Schematic representation of (a) the classical SK 
growth, and (b) the coherent SK growth. In the latter case the 
side walls are shown steeper to demonstrate the compression 
exerted by the wetting layer. The MDs in (a) are denoted by 
inverse T's. 

FIG. 2. The pairwise potential of Eq. (1) with /i — 12, 
v = 4 and Vo ~ 1- The dashed vertical line through the 
inflection point Xi separates the regions of distortion {x > Xi) 
and undistortion {x < Xi) of the chemical bonds shown in the 
upper part of the figure. 

FIG. 3. Schematic view of multilayer islands with different 
slopes of the side walls (a) 60°, (b) 30°, and (c) 19.1°. 
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FIG. 4. Illustration of the solutions of the one-dimensional 
model of Frank and van der Merwo: (a) a chain without a mis- 
fit dislocation, (b) a misfit dislocation in a compressed chain, 
(c) a misfit dislocation in an expanded chain. With increasing 
the chain length in (a) the end atoms are more displaced from 
the bottoms of the potential troughs and approach the crests 
between them. 

FIG. 5. Distribution (a) of the strain — ^i+i — — f 
in a monolayer height compressed (/ = 0.07) and expanded 
(/ = —0.07) chains, and (b) of the corresponding bond energy 
in units of K. W/Vo = 1/3, ^ = 12, i/ = 6. 

FIG. 6. Distribution of the adhesion energy 'Pi/Vo — 1 in 
monolayer height compressed (curve 1) and expanded (curve 
2) chains. W/Vo = 1/3, ^=\2,v = 6. 

FIG. 7. The mean adhesion parameter # as a function 
of the number of atoms in the chains for positive (/ = 0.07) 
and negative (/ = —0.07) values of the misfit. W/Vo = 1/3, 
1^ = 12, iy = 6. 

FIG. 8. Distribution of the total energy (strain plus adhe- 
sion) in units of Vo in monolayer height compressed (/ = 0.07) 
and expanded (/ = -0.07) chains. W/Vo = 1/3, // = 12, 
v = 6. 

FIG. 9. Energy per atom of pyramidal 3D islands in units 
of Vo with different slopes of the side walls denoted by the 
figures at each curve as a function of their thickness in number 
of monolayers. The number of atoms A'^i = 19 in the base 
chain is one and the same for all curves. The frustum of the 
pyramid with a slope of 60° of the side walls and height of 9 
monolayers represents the equilibrium shape. W/Vo = 1/3, 
^ = 12, = 6. 

FIG. 10. The dependence of the energy per atom on the 
total number of atoms in compressed coherently strained is- 
lands with different thickness in monolayers denoted by the 

figures at each curve: (a) / — 0.03, (b) / = 0.05. The curve 
denoted by MD in (b) represents the energy of a monolayer 
chain containing one misfit dislocation. The numbers 

N23, etc. give the limits of stability of monolayer, bilayer 
islands, respectively. W/Vo ~ 1/3, = 12, = 6. 

FIG. 11. The dependence of the energy per atom on 
the total number of atoms in units of Vo in expanded coher- 
ently strained islands with different thickness in monolayers 
de- noted by the figures at each curve, and at large negative 
value of the misfit / — —0.1. The curve denoted by MD rep- 
resents the energy of a monolayer chain containing one misfit 
dislocation. W/Vo = 1/3, n = 12, u = 6. 

FIG. 12. Misfit dependence of the critical size Ni2- The 
vertical dashed lines denote the critical misfits below which 
N12 is infinite. The curves are shown in one quadrant for 
easier comparison. 



FIG. 13. The energy change AS„ in units of Vo connected 
with the transformation of mono- to bilayer islands (curve 

1-2), bi- to trilayer islands (curve 2-3), and of tri- to fourlayer 
islands (curve 3-4), respectively, as a function of the number 
of atoms n in the uppermost chain. / = 0.05, W/Vo = 1/3, 
^ = 12, = 6. 
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